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Abstract. A model in which a gap forms in the renormalized electronic density of 
state (DOS) with missing states recovered just above the pseudogap A p9 , is able to 
give a robust description of the striking, triangular like, peak seen in the real part 
of the optical self-energy of underdoped cuprates. We use this model to explore the 
effect of the pseudogap on the real part of the optical conductivity and on the partial 
sum rule associated with it. An important result is that the optical spectral weight 
redistributes over a much larger frequency window than it does in the DOS. 
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1. Introduction 

Cuprate superconductors undergo considerable change in their electronic structure as 
a function of doping. In the overdoped regime there is evidence that Fermi liquid 
theory applies. However, as the doping is reduced through the optimum and then 
underdoped regime, a non-Fermi liquid state evolves which is characterized by the 

formation of a pseudopgap[H ^SHIQIQHSISHniQIlQS. There is as y et no 
consensus as to the correct microscopic understanding of the pseudogap state. There is 
evidence that it corresponds to the formation of preformed Cooper pairs 0] at some new 
characteristic temperature T* above the superconducting T c . Superconductivity sets in 
at the temperature where phase coherence among the preexisting pairs occurs. Only 
then does the system acquire long range order. Another prominent possibility for the 
pseudogap is competing order [III [13 IIEl IIZl 03 OS EIB Ell E21 1231 Ell 123 ESI EZl I2B] 
such as d-density wave formation. On the other hand, from a phenomenological point of 
view many experimental findings, in the pseudogap region of the cuprate phase diagram, 
can be understood in a model which corresponds to a reduction of the electronic density 
of state in the vicinity of the Fermi energy over an energy region (A pg ) which sets the 
pseudogap scale. An important concept in such a characterization of the pseudogap is 
its variation with temperature. At zero temperature there is a full gap of order A pg 
which fills in with increasing temperature but does not change its magnitude^ [12]. 
The pseudogap temperature (T*) corresponds to compete filling rather than closing i.e. 
A pg - 0. 

Recently more details about the temperature evolution of the pseudogap have 
emerged from angle-resolved photoemission (ARPES) experiments which have been 
interpreted in term of Fermi arcs[7j. Below T* a full pseudogap opens up, but only 
on a small region of the Fermi surface near the antinodal direction. The rest of the 
Fermi surface, called the Fermi arc, remains ungaped. Experiments have shown that 
the length of the Fermi arc centered in the nodal direction is proportional to the reduced 
temperature t = T/T* and vanishes at T = at which point the entire Fermi surface is 
fully gaped. The existence of a full gap in the electronic density of state at T = 
is consistent with many other experiments in particular specific heat. It also has 
implications for the behavior of the in-plane optical conductivity. 

It has recently been pointed [29J out that the sharp triangular like cap observed in 
the real part of the optical self-energy seen in underdoped samples of Bi 2 Sr 2 CaCu208+<5 
(Bi-2212)[Hl EI] and ortholl YBa 2 Cu 3 6 . 5 (YBC0 6 . 5 )l32j follows directly from the 
opening of a full pseudogap with the lost electronic density of states below A pg 
recovered in the region just above it. This represents a clear signature of pseudogap 
behavior in optical spectroscopy. In this paper we consider the implication of such a 
phenomenological model on optical properties. 

The paper is structured as follows. In section II we introduce the generalized Drude 
model for the optical conductivity which relates it to an optical self-energy, E op (u;). We 
also summarize the data for the real part of E op (a>) on which our pseudogap model 
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is based and the distinct difference between underdoped and overdoped cases is noted 
and emphasized. Our theoretical model is introduced and compared with the data. In 
section III we deal with signatures of the pseudogap in the real part of the conductivity 
and partial sum rule. In section IV we conclude the paper with a summary of our 
findings. 



2. Theoretical model 



The optical conductivity c(uj) in a correlated electron system can be analyzed in terms 
of a generalized Drude form written as [291 [30] 



n 2 

" (T '" ) = ^- 2 ^(r>r (1) 

where T is temperature, Q p is the plasma frequency and £ op (T, uj) = S° P (T, uj) + 
iYf^iT, uj) is the optical self-energy. The imaginary part of -2£ op (T, uj) is equal to the 
optical scattering rate l/r op (T, uj) and the real part can be written in terms of an optical 
effective mass m*' op (T, u)/m with u[m*' op (T, uj)/m—l] = — 2E° P (T, uj). While the optical 
scattering rate and the mass renormalization X op (T, uj) (1 + X op (T, uj) = m*' op (T,uj)/m) 
defined here are not the same as those defined from the quasiparticle self-energy 
£ 9P (T, uj) they are related through the equation for the conductivity. Neglecting vertex 
corrections and taking zero temperature (T = 0) for an isotropic system we have [3 3 j 

SI 2 i r u i 
a ^~4^uJ dU uj + i/r imp - }Ji p (uj') - ^ p (u - u') 

where we have also included the possibility of elastic impurity scattering through the 
constant scattering rate 1/T imp . The quasiparticle scattering rate is — 2J2f(uj) and 
uj[m*' qp (uj)/m — 1] = — 2£f (u;) in complete analogy with the optical case. If a boson 
exchange theory is used to describe the interactions among electrons, the quasiparticle 
self-energy at T = is related to the electron-boson spectral density I 2 x{^) through 
the equation [M], |35] 



r°° q — uj 

E qp (uj)= dttl 2 x (ty\n ml dttl 2 x (ty (3) 

Jo il + uj J Q 

Here we need to generalize the formalism just given to include the possibility of an 
energy dependent renormalized electron density of state (DOS) N(u) which is defined 
as [351 E6J S3 EH [39] 

N{ U ) = y- ImG ^ u) (4) 
k 

where G(k, uj) is the fully renormalized Green's function. In this case Eq. [3] needs 
modification and reads [31 [351 M, EZJ [38j [39] 

Z qp (u)=2uP duj'N(uj') dtt— ^ — in / dttI 2 X (n)N(uj-tt)(5) 

Jo Jo uJ 2 -(uj' + iiy J 




Figure 1. (color online) Minus twice the real part of the optical self-energy — 2E° p (w) 
in units of cm -1 as a function of oj also in cm" 1 for an underdoped Bi-2212 with T c = 
69 K (UD69) and an overdoped Bi-2212, T c = 60 K (OD60) at two temperatures; 
dashed (blue) in superconducting state and sold (red) in normal state. Shown in the 
inset are equivalent results for ortholl YBCO6.50 which is underdoped with T c = 59 K 
(UD59). 



Here N(u) = [N(u) + N(—u)]/2 is the symmetrized DOS. For the case with an energy 
dependent DOS the relationship between S 9P and S op is more complicated than Eq. 
EJ There exists however simplified equation for X° p (uj) and 1/t° p (uj) which, while not 
exact, are sufficiently accurate for the present discussion. They are[ 



2 f°° f°° r (t / 4- O") 2 i 



and 

1 2tt 



r2„ 



. ^/ 2 X(^) / du/iV(aO (7) 

t°p(^) w Jo io 

In this approximation, optical and quasiparticle quantities are related by YP p (uj) = 
d / duj[u)Y> op (u))] which can be verified through direct differentiation of Eq. [6] and and 
comparison with Eq. [5j This relationship has been used in a recent comparison of high 
energy scales seen in optical data with those seen in ARPES[l0] which measures directly 
the quasiparticle self-energy. 

The exact microscopic origin of the pseudogap is not known. Here we model it as 
a gap in the fully renormalized electronic density of state N(u) of formula Eq. HI Such 
a model has been used previously to analyze the specific heat[10l[TT] in the underdoped 
regime of the high T c cuprates and more recently applied to optics[29]. The motivating 
data is reproduced in Fig. [T]for the convenience of the reader. These sets of data, for the 
real part of the optical self-energy are presented for an in-plane underdoped sample of 
Bi-2212 with a T c of 69 K, for another overdoped sample with T c = 60 K and in the inset 
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Figure 2. (color online) Model calculation of minus twice the real and imaginary parts 
of the optical self-energy — 2S^ p (w) (top frame) and — 2Y,2 P (u) (bottom frame) in cm -1 
as a function of u> also in cm . All curves are based in the electron-boson spectral 
density shown in the inset lower frame. The dashed (blue) curve has no pseudogap 
while the dash-dotted has a pseudogap A pg — 550 cm" 1 but with no recovery of states 
in the self consistent density of state N(w). The others have a recovery region right 
above A pg with conservation of total states applied. The recovered states are piled up 
in the region A p9 , (1/2) A pg , and (1/10) A M for medium (purple), light (olive), and 
heavy (red), respectively. The inset in the upper frame shows the effective DOS N(cu) 
in the case (1/2)A PS . 

data on underdoped ortholl YBCO6.50 with T c = 59 K. This material is particularly well- 
ordered with every second chain full and the others completely empty. In all cases two 
values of temperature are shown, one in the superconducting state (dashed blue curve) 
and the other in the normal state just above T c (solid red curve). The difference between 
the behavior of underdoped and overdoped samples is striking and can be understood [29J 
as due to the opening of a pseudogap in the fully renormalized density of state N(u) 
of formula Eq. H] which also determines the optical self-energy through equation Eq. [6] 
and Eq. [71 The prominent peak around 750 cm -1 seen in both underdoped materials 
which is absent in the overdoped case, can be traced to the opening of a gap in N(u) 
with lost states piled up in the energy region just above u = A pg as well as the existence 
of a prominent boson mode in the electron-boson spectral density I 2 x{ty- All three of 
the above conditions are needed. 

Fig. [2] shows results of model calculations for minus twice the real part (top frame) 
and imaginary (bottom frame) part of the optical self-energy (E op (u;)) for a system 
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with the electron-boson spectral density I 2 x{oj) shown in the inset of the lower panel. 
This model for the spectral density (note the prominent resonance peak at 31 meV) is 
motivated by an earlier study [32] of ortholl YBCO6.50 in which I 2 x{oj) was fit to its 
measured optical scattering rate. What is different here is that we includes a pseudogap 
in N(uj) as well as a recovery region just above A pg where the lost states in N(uj) are to 
be found so as to conserve states in the DOS. The blue dashed curve is for comparison 
and includes no pseudogap. The expressions for \° p (uj) and \/t op {uj) in this simple limit 
reduce to [31] 

\°p(uj) = - / dtti 2 x (n) 

w Jo 

and[Ml SB 12] 

— — = - / dn(u - n)i 2 x ({i) (9) 

T op (uj) UJ J 

On comparing Eq. with the imaginary part in Eq. [31 it is clear that for a simple 
delta function 7 2 x(f2) = I 5(Q — Qe) which represents coupling to a single Einstein mode 
at Q E , the optical scattering rate starts out of zero at uj = Q E , then rises according 
to the factor (uj — D,e)/oj and reaches its saturated value of 2nl only for uj 3> fl E - 
By contrast the quasiparticle scattering rate has a discontinuous jump out of zero at 
uj = Q E to its saturated value and remains at this constant value for all energies beyond 
this. This behavior for 1/t° p (uj) is seen in the (blue) dashed curve of Fig. [2] bottom 
frame. Because we have also included a background in I 2 x{^) in addition to a prominent 
peak at Qe = 250 cm -1 , there are minor differences, including very small tails below the 
energy of the prominent peak in I 2 x{^>) shown in the inset. The dot-dashed (blue) curve 
includes a pseudogap A pg = 550 cm -1 with the lost states in N(uj) moved to infinity. In 
this case the main rise in 1/t° p (uj) is at Qe + A pg ~ 800 cm -1 . Beyond this the curve 
rises approximately like [uj — {VLe + A. pg )]/uj which is less rapid than the (uj — VLe)/uj 
curve for the A pg = case. 

The three remaining curves also have a pseudogap of 550 cm -1 but, in addition the 
lost state in N(uj) are placed in the energy region just above A pg and state conservation 
is respected. This causes the scattering rate 1/t° p (uj) to rise much faster than in the 
dot-dashed curve because just above the gap we have more states to scatter into. The 
steepness of the rise depends on the distribution of states above uj = A pg . The more 
the pileup is restricted in range the steeper the rise. Medium (purple), light (olive), and 
heavy (red) solid curves correspond, respectively to the case when the missing states 
are placed between A pg and 2A p9 , 1.5A ra , and l.lA pg (see inset in the top frame for 
the case 1.5 A pg where the DOS is shown). Note also that the end of the recovery 
region in all cases is mark with a kink in 1/t° p (uj) after which the remaining rise is 
much more gradual and smooth. The features just described imply definite signatures 
in the corresponding real part of the optical self-energy as these are related by Kamers- 
Kronig (K-K) transform. The results for — 2E° p (o;) based on Eq. 6 are shown in the top 
frame of Fig. [21 As is known from the work of Carbotte, Schachinger and Hwang [34J 



^In 

UJ 



UJ 



O 2 



+ In 



Q — uj 



(8) 
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the dashed (blue) curve would have a logarithmic singularity in slope aX u = VL E ii 
we were using a pure delta function model and would have zero slope at uj = \^2Qe- 
Similarly the dash-dotted (blue) curve would have infinite slope at VLe + A pg and zero 
slope at v2(Qe + A pg ). These rules are very nearly satisfied in our model calculation 
even though we are using the spectra displayed in the inset lower frame rather than a 
pure delta function. We recall that the K-K transform of a sharp step like rise at Qe 
and constant after this as applies to the quasiparticle scattering rate in a delta function 
model, has a logarithmic singularity at uj = Qe- The heavy (red) curve in the lower 
frame comes close to this case and indeed its K-K transform shows a sharp peak at this 
frequency reminiscent of a logarithmic singularity. We believe this to be the signature in 
the real part of the optical self-energy of pseudogap formation [29J as seen so prominently 
in the data of Fig. [T] for the two underdoped samples. A model with recovered DOS 
within A pg shows a clearly identifiable hat type structure in the medium solid (purple) 
curve for — 2S° p (u;) (Fig. [2]top frame) missing in both (blue) curves. This hat is perhaps 
not quite as sharp in this model calculation as it is in the data which is however less 
peaked than the heavy (red) curve. This indicates that the recovery region is consistent 
with a renormalized density of state for which the conservation of states occurs on the 
scale of < A pg . 

Another interesting possibility to consider is the case when the pseudogap does not 
reduce N(uj) to zero for u < A pg but rather still has a finite value. To illustrate this 
possibility, in the top left hand frame of Fig. 3, we show the imaginary part of E op (u;) 
for the case when the DOS is reduced to half its value rather than to zero below A pg . 
The light solid (olive) curve is to be compared with the heavy solid (red) curve which 
we reproduced from the bottom frame of Fig. 2. Both curves include full recovery in the 
energy interval A pg and 2A P9 . The light solid (olive) curve now starts at uj = Qe and 
has a step at Qe + A p9 after which it shows the characteristic sharp rise which we have 
associated with the density of state recovery region. In this case we see clear signatures 
of the resonance mode and of the pseudogap recovery region separately. 



3. Effect of the pseudogap on real part of conductivity and on its partial 
sum rule 

Having a model for the optical self-energy H op (uj), we can calculate from it the real part 
of the conductivity. From Eq. [Hand the definitions of 1/t° p (uj) and X° p (uj) we get 

a (uj) = ^ p 1/T ° PM (10) 

XK ' Air [uj(1 + \°p(uj))} 2 + {1/t°p(uj)] 2 1 ' 

One point that needs to be made is that at zero temperature the inelastic optical 
scattering rate vanishes for uj < Qe in an Einstein model and the optical conductivity 
Eq. [10] becomes pathological: Oi(a>) = Q p 5(uj)/[A(1 + A op (0))] has a delta function at 
uj = with weight reduced over the free electron case by 1/(1 + A), where A = A op (0), 
which coincides with the quasiparticle mass enhancement parameter X qp (uj = 0) even 
for the case of an energy dependent density of state as we are considering here. . This 
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can be remedied by including a small amount of elastic impurity scattering as we did 
in Eq. [2l In this case 1/t° p (uj) around u = is a finite constant and its corresponding 
real part vanishes. Results for o~i(u) vs. uj are shown in the bottom frame of Fig. [3j In 
both cases \/r imp is set equal to 80 cm -1 and Q p is 10000 cm -1 . For the left hand frame 
we have used the I 2 x{^) shown in the inset of the bottom frame of Fig. [2] which has 
an optical resonance as well as a background while for the right hand frame only the 
background is used. Starting with the left panel, the dashed (blue) curve clearly shows 
two regions, the Drude plus a Holstein boson assisted absorption piece which extends 
way beyond the Drude and mirrors the spectral density I 2 x{^)- For a delta function 
its onset is at uj = Qe arid grows out of zero according to Eq. [TUl It contains A/(l + A) 
of the optical spectral weight. The remaining optical spectral weight 1/(1 + A) is to be 
found in the coherent Drude contribution. The width of the Drude is given very nearly 
by [l/r imp (l + A)] which is different for the various curves because A varies as we will 
explain. In the inset (bottom left frame) we repeat the solid (red) curve for pseudogap 
plus recovery and compare it with the same case but now the elastic scattering has been 
increased by a factor of 4. This fills in the region between Drude and Holstein processes 
much as is seen in the Bi-2212 data as one goes further into the underdoped region to 
T c = 69 K[3T]. This can effectively switch some of the extra spectral weight that was 
transferred to the coherent part of the conductivity (Drude part) by the opening of the 
pseudogap back to the larger energies associated with the incoherent boson assisted part 
of the absorption (Holstein) region. 

When there is no pseudogap as in the dashed (blue) curve, the quasiparticle mass 
enhancement parameter is given by 



If however a pseudogap is introduced it becomes modified and reads instead from Eq. 
[5] in the limit uj — ► 



beyond. That is, we have piled up the missing state equally in the interval A pg to 2A, 
In this case the mass enhancement parameter becomes 




(11) 




(12) 



To remain simple we assume N(u) = t for uj < A pg and 2 — t for uj G (A pg , 2A pg ) and 1 




(13) 




(14) 



(15) 
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Figure 3. (color online) Model calculation of minus twice the imaginary part of 
the optical self-energy — Ttf^iio) (top frame), the partial optical sum rule to energy 
u> in units of Sl^cm -2 (middle frame) and the real part (absorptive) of the optical 
conductivity in cm -1 (bottom frame) as a function of w in cm -1 . The curves on the left 
were calculated with the I 2 x(^>) shown in the inset of Fig. 2 while the right hand panel 
employed only the background contribution without a sharp optical resonance mode. 
Dashed (blue) curve is for no gap, dash-dotted (purple) for a gap without recovery, 
solid (red) for a gap with recovery, and light solid (olive) for a model in which only 
half the states are removed below the gap (with recovery). Inset same as solid (red) 
curve but with 4 times more elastic (impurity) scattering included. 



for the case of no recovery region. This factor in effect reduces the low Q contribution to 
the mass enhancement factor from the spectral density alone. For a full pseudogap the 
suppression provides an extra factor of 3Q/(2A pg ) and Q/A pg at small Q while for finite 
t it gives t + 3Q/(2A pg ) and t + Q/A pg , respectively For the specific case considered 
here A = 2.37, when we include a pseudogap Eq. [12]is reduced to 0.78 for a full gap with 
no recovery region above it and it is 1.06 when recovery is included as in Eq. [TH An 
unexpected consequence of this reduction in A is that the Drude contribution (coherent 
part) to <J\{uS) is increased when a pseudogap is included and there is a corresponding 
decrease in the boson assisted Holstein (incoherent) contribution centered off u — 0. 
These factors translate into a wider Drude as A is decreased. Dash-dotted (purple) 
is widest, then solid (red) and finally dashed (blue) at the same time the Holstein 
boson assisted region shows increasing weight from dash-dotted (purple) to solid (red) 
to dashed (blue) curves. Note also that a complete pseudogap below uj = A pg cuts 
off the Holstein region which now starts at u = A pg + Qe in the both solid (red) and 
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dash-dotted (purple) curves. The light solid (olive) curve is different from the others 
and shows two peaks rather than one in the Holstein region. It corresponds to an 
incomplete pseudogap with height t = 0.5 at uj = 0. This translates into a larger A 
than for the solid (red) curve and its Drude peak indeed falls between (red) solid and 
(blue) dashed curves. Its Holstein region however starts at Qe — 250 cm -1 (position of 
large peak in J 2 x(u;) of the inset in lower frame of Fig. [2]) because the DOS is finite at 
the Fermi surface. This onset is followed by a peak with a second onset seen clearly at 
uj = A pg + Qe which corresponds to the sudden increase in the electronic DOS. To end 
this discussion of the effect of a pseudogap on the real (absorptive) part of the optical 
conductivity we consider the lower right hand panel of Fig. 3. Here only the background 
spectrum in I 2 x(oj) is included i.e. the delta function like optical resonance is excluded. 
We see that now Drude and Holstein contribution in the dashed (blue) curve are not 
clearly separated. While a two-contribution structure is seen in the other two curves 
these are not as well defined as in the corresponding curves of the left-hand panel. The 
existence of a sharp peak in electron-boson spectral density helps separate out the two 
distinct absorption processes (Drude and Holstein). 

In the middle frame of Fig. 3 we show results for the partial optical sum S(uj) 
defined as 

S{u) = / duj'a^uj'). (16) 
Jo 

For uj — > oo in Eq. [16j we get the usual complete sum rule with S(uj — > oo) = fip/8. 
In our units this corresponds to ~ 2.6 x 10 6 cm~ 2 with Q p = 10000 cm -1 . One sees 
clearly in these curves a rapidly increasing Drude contribution followed by a flattened 
region and then the Holstein contribution setting in at large values of uj. The fraction 
of spectral weight seen in the Drude region is 1/(1 + A) of the total contribution. The 
remainder A/(l + A) is the boson assisted contribution. If these two regions were truly 
separated in cf\{uj) the flattened region noted above would be perfectly flat. We make 
one more point about these results. In all cases there is, of course, conservation of 
optical spectral weight but this conservation occurs over a much larger frequency region 
for the real part of the conductivity than for the density of states itself which, in our 
model, is limited to the range to 2A pg . This is also true of the imaginary part of the 
optical self-energy but not for ct^uj) or the partial sum S(u) of Eq. [T5j Yu et a/.|i3] 
have recently noted this in their c-axis optical study of spectral weight redistribution 
due to the pseudogap. 

Finally returning to the top frame of Fig. 3 we note that the sharp onset in 
scattering rate seen in the left-hand panel is considerably smeared out when the resonant 
peak in the boson spectral density I 2 x{^) shown in the inset of the lower frame of Fig. 
2 is left out and only the background spectral density is employed. This is shown in the 
right-hand panel (compare color coded curves). 
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4. Conclusions 

Motivated by the observation of hat like peak structures seen in the real part of the 
optical self-energy in underdoped cuprates which are absent in the overdoped case, we 
have considered a pseudogap model for their electronic structure consisting of a simple 
gap (A pg ) in the electronic DOS about the Fermi energy, with missing states recovered 
in the energy region right above it. This simple model augmented with an optical peak 
in the electron-boson spectral density is remarkably successful in describing the data 
and is further used to describe other optical quantities. The optical scattering rate, 
the real part of the conductivity and its partial sum are considered as is the derivative 
of E° p (c<;). These quantities all show specific signatures of the pseudogap. While the 
redistribution of electronic states in the DOS is limited to the region ui < 2A p9 we found 
that the corresponding redistribution of optical spectral weight in the real part of the 
conductivity is spread over a much large range. This is also true for the real part of the 
optical self-energy but not for its imaginary part for which the important changes due 
directly to the opening of pseudogap are confined much more to the range u < 2A pg . 
Comparing results of model calculations we conclude that pseudogap effects can be 
readily distinguished from the effect of multiple peaks in the electron-boson spectral 
density 7 2 x(cj). 
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